Orbital entanglement and violation of Bell inequalities in mesoscopic conductors 
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We propose a spin-independent scheme to generate and detect two-particle entanglement in a 
mesoscopic normal-superconductor system. A superconductor, weakly coupled to the normal con- 
ductor, generates an orbitally entangled state by injecting pairs of electrons into different leads of 
the normal conductor. The entanglement is detected via violation of a Bell inequality, formulated 
in terms of zero-frequency current cross-correlators. It is shown that the Bell inequality can be 
violated for arbitrary strong dephasing in the normal conductor. 

PACS numbers: 74.50.+r, 73.50.Td, 73.23.Ad 



Entanglement is one of the most intriguing features 
predicted by quantum theory It leads to correlation 
between distant particles, which can not be described by 
any local, realistic theory This nonlocal property 
of entanglement has been demonstrated convincingly in 
optics where entangled pairs of photons have been 
studied over several decades. Apart from the fundamen- 
tal aspects, there is a growing interest in using the prop- 
erties of entangled particles for applied purposes, such 
as quantum cryptography and quantum computation 

Recently, much interest has been shown for entangle- 
ment of electrons in solid state devices. A controlled 
generation and manipulation of electronic entanglement 
is of importance for a large scale implementation of quan- 
tum information and computation schemes. Electrons 
are however, in contrast to photons, massive and elec- 
trically charged particles, which raises new fundamental 
questions and new experimental challenges. Existing sug- 
gestions are based on creating jfj, 0] , manipulating and 
detecting H,EJEi spin-entangled pairs of electrons. This 
requires experimental control of individual spins via spin 
filters or locally directed magnetic fields on a mesoscopic 
scale. Here we propose a spin-independent scheme for 
creating and detecting orbital entanglement in a meso- 
scopic normal-superconductor system. 

We show that a superconductor, weakly coupled to a 
normal conductor (see Fig. H), creates an orbitally en- 
tangled state by emitting a coherent superposition of 
pairs of electrons into different leads of the normal con- 
ductor. The zero-frequency correlation between currents 
flowing into different normal reservoirs is shown to be 
equivalent to a pair coincidence measurement: only cor- 
relations between the electrons from the same entangled 
pair contribute. As a consequence, a standard Bell In- 
equality (BI) can be directly formulated in terms of the 
zero- frequency current correlators. We find that a vio- 
lation of the BI, demonstrating the entanglement of the 
pair state, can be obtained for arbitrary dephasing in the 
normal conductor. 

We first consider a simplified version of the system 
(see Fig. HJ, a more detailed discussion is given below. 



A single superconductor (S) is weakly coupled to 
a normal conductor, a ballistic two-dimensional electron 
gas, via two tunnel barriers 1 and 2 with transparency 
r < 1. The normal conductor consists of four arms, 
1A, IB, 2A and 2B, with equal lengths L. The arms 1A 
and 2A (IB and 2B) are crossed in a controllable beam 
splitter A(B), parameterized via the angle ^a{.^b), an d 
then connected to normal reservoirs +A and —A (+B 
and — B). The beam splitters ^3] are assumed to support 




FIG. 1: The system: A single superconductor (S) is connected 
to four normal arms via two tunnel barriers 1 and 2 (thick 
black lines). The arms are joined pairwise in beam splitters 
A and B and end in normal reservoirs + and — . 

only one propagating mode. The states 1+,??) and | — ,rj) 
for electrons going out into the normal reservoirs and 
the states |l,ry) and \2,rj) of the electrons emitted from 
superconductor are related via a scattering matrix: 
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, V = A,B.(1) 



The angles 4>a and <\>b can be tuned between and 7r/2 by 
tuning the beam splitter from "open" , when the electrons 
are transmitted through from 1(2) to — (+), to "closed", 
when the electrons are completely reflected, from 1(2) to 
+ (— ). We consider the low temperature limit, kT <C eV. 
A negative voltage — eV is applied to all the normal reser- 
voirs and the superconductor is grounded. The voltage 
eV is smaller than the superconducting gap A, so no 
single particle transport takes place. It is further as- 
sumed that the size of the system is smaller than the 
phase breaking length. 
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We first present a simple and transparent explanation 
of how the entanglement is generated and detected, a 
rigorous derivation follows below. The superconductor 
emits pairs of particles into the normal arms. Since the 
superconductor is a single object, coherent on a macro- 
scopic scale, the state of an emitted pair is a linear super- 
position of states corresponding to a pair emitted through 
barrier 1 and 2. The emitted pair can either split, with 
one electron going to each beam splitter A and B, or both 
electrons can go to the same beam splitter. However, the 
latter process does not contribute |l3( to the quantity of 
interest, the current cross-correlations, to leading order 
in pair emission (Andreev reflection) probability, propor- 
tional to r 2 . The relevant part of the state of the emitted 
pair can thus be written as 



1*15 



I* 



AB 



1*12 



(|11) + |22))/V2, (2) 



a product of a state 1*12), orbitally entangled with re- 
spect to emission across barrier 1 and 2 (called 12-space 
below), and a state |*ab)> describing one electron going 
towards A and one towards B, containing all additional 
information, such as energy and spin dependence. The 
beam splitters A and B, parameterized via the angles 4>A 
and <f>B, rotate the state 1*12) [see Eq. JD]. 

The entanglement is detected by violation of a BI. We 
recall that in the original formulation Q , a source emit- 
ting spin-1/2 sing lets was considered. The BI, as formu- 
lated in Ref. 1141 . 
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is expressed in terms of spin correlation functions [if 



E((j> A ,c]> B )=P+ 



P+--P-+ + P— 



(4) 



Here P a /3{4>A, 4>b) are the joint probabilities to observe 
one particle in detector A with a spin a = ± [+(— ) de- 
noting up (down)] along the ^-direction, and the other 
in detector B with a spin (3 = ± along the 0s-direction. 
The joint probabilities are given by 



Pa/3 {4>A,<t>B) = 1 + a/3 COs[2((pA ~ 4>b)}- 



(5) 



Inserting the values for P a ^ from Eq. (JfTJ into Eq. J3Jl, 

we get E{(j>At<l>B) = cos[2(<^4 — <j) B )\. We then directly 
see that for e.g. angles (f>A — ir/8,4> B — tt/4, (j)' A = 37r/8 
and 4>' B = tt/2, the BI in Eq. © is maximally violated, 
i.e. we get S = 2\f2. 

In our orbital setup (see Fig ^) , it is clear from Eq. 
(0) that the 12-space plays the role of a pseudo-spin 
space and the normal reservoirs act as detectors. We 
can thus, in direct analog to Ref. 2], formulate a BI in 
terms of any observable which is directly proportional to 
the corresponding joint probability P a f3{4>A, <Pb) for our 
state 1*12) [here a,/3 = ± denote the reservoirs, see Fig. 



We find below that the zero- frequency current cross- 
correlator is given by 

/oo 
dt{8i aA {t)8i 0B {Q)) = PqPcc^aAb), (6) 
-00 

i.e. directly proportional to the joint probability distribu- 
tion. Here P = e 3 Vr 2 /2h, and 6I an (t) = I av {t) - (I ari ) 
is the fluctuating part of the current I ar i(t) in reser- 
voir arj. This leads to the important result that the 
Bell inequality in Eq. (0 can be directly formulated in 
terms of the zero-frequency current correlators in Eq. 
©■ We note 0] that when inserting S a p directly into 
the correlation functions Eq. J3} (i.e. substituting P a p 
with S a p), we must divide by the sum of all correlators 
S ++ + +S |_ + <!?__ = Pq, which just eliminates Pq. 

The simple result in Eq. © can be understood by 
considering the properties of the time-dependent corre- 
lator (S I a A {t)5I[]B (0)). It is finite only for times t < r c , 
where r c = Ti/eV is the correlation time of the emitted 
pair (see Fig. |5J). In the tunneling limit under consid- 
eration, r < 1, the correlation time is much smaller 
than the average time between the arrival of two pairs 
e/I ~ h/eVT 2 . As a result, only the two electrons within 
a pair are correlated with each other, while electrons in 
different pairs are completely uncorrelated. Thus, the 
zero frequency current correlator in Eq. © is just a coin- 
cidence counting measurement running over a long time, 
collecting statistics over a large number of pairs [R| . 

For a rigorous derivation of the above result, we first 
discuss the role of the superconductor as an emitter of 
pairs of orbitally entangled electrons. Since the system 
is phase coherent, we can work within the scattering ap- 
proach to normal-superconducting systems 
starting point is the many-body state of the normal reser- 
voirs, describing injection of hole quasiparticles at ener- 
gies from to eV (suppressing spin notation), 

1*-)= II lUEhl A (Ehl B (Ehl B (E)\0), (7) 

Q<E<eV 

where the ground state |0) is the vacuum for quasiparti- 
cles in the normal reservoir. The operator j\ A (E) creates 
a hole plane wave at energy E (counted from the super- 
conducting chemical potential /zs) in lead 1A, going out 
from the normal reservoirs towards the superconductor, 
and similarly for the other operators. The commutation 
relations for the 7-operators are standard fermionic. 

To obtain the state of the quasiparticles going out from 
the superconductor, |* D ut), we note that the operators 
creating and destroying quasiparticles going out from the 
superconductor (towards the normal reservoirs) are re- 
lated 0, 0| to the operators of the incoming quasipar- 
ticles via a scattering matrix. The scattering amplitude 
for a hole injected in arm 1A, to be back-reflected as an 
electron in arm 2A, is denoted 1A and similarly for 
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the other scattering amplitudes. In the tunnel limit un- 
der consideration, the amplitude to backscatter as the 
same type of quasiparticle is ^™ ., , ~ 1- The amplitude 
for Andreev reflection, r|^-/„/, is given by — ir/4 (in- 
dependent on energy). Processes where a hole incident 
on barrier 1(2) is backscattered as an electron at barrier 
2(1), i.e. when a pair in the superconductor breaks up, 
are exponentially suppressed with the distance between 
the two emission points [|J and can be neglected in the 
present setup. 

The tunneling limit T <C 1 makes it relevant to change 
the perspective from a quasiparticle picture to an all- 
electron picture. An Andreev reflection, occuring with 
a small probability T 2 , can be considered as a pertur- 
bation with respect to the ground state |0) in the nor- 
mal reservoirs (a filled Fermi sea of electrons at energies 
E < — eV). It creates an excitation consisting of a pair of 
electrons (see Fig. - Formally, performing a Bogoliubov 
transformation (i.e. 7^ (E) = c] v (-E), with c] v (E) be- 
ing a standard electron creation operator), the state of 
the quasiparticles going out from the superconductor be- 
comes |2Jj] to first order in Y 



|* out > = |0> + |*> + |vl/), 
where the states 



(8) 
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c\ A (E)c\ B (-E) 
+ cl A (E)4 B (-E)] |0), 



(9) 



describe orbitally entangled electron "wave packet" pairs, 
i.e. a superposition of pairs at different energies [2lj |. 
In a first quantized notation, the state \^) is just the 
state in Eq. (0). We emphasize that the change from a 
quasiparticle to an all-electron picture, providing a clear 
picture of the entanglement, does not alter the physics. 

The detection of the entanglement is done via the 
zero frequency current cross-correlators S a p in Eq. 
®. To obtain 

S a p j we insert the sec ond quantiza- 
tion expression for the current operator [l9j |. I av (t) = 
(e/h) J dEdE'exp[i(E - E'jt/h^iE'jc^iE), into Eq. 
© and average with respect to the state \^) in Eqs. 
JHJ and @. The beam splitters are taken into ac- 
count by relating the c^-operators in arms 1A and 2A 
(IB and 2B) to the c Q?) -operators in the reservoirs +A 
and —A (+B and — B) via the scattering matrix in 
Eq. The average current, equal in all arms arj, 

is (I av ) = I = (e 2 /2h)T 2 V , independent of the beam 
splitter transparency. 

The ground state |0) in Eq. JHJ does not contribute to 
the correlator. Moreover, calculating the current corre- 
lator we find that the state | ^) , describing two electrons 
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FIG. 2: Left, upper: In the filled stream of incoming holes 
(empty circles) from the normal reservoir, occasionally a hole 
is back-reflected as an electron (solid circles). The "miss- 
ing hole" (i.e an electron) and the Andreev reflected electron 
constitute the pair emitted from the superconductor. Left, 
lower: The correlation time r c = fi/eV (width of the wave 
packet) and the average time between emission of two subse- 
quent pairs e/I = Ti/eVT 2 3> t c of the current correlator. The 
small time difference 7i/A <C t c between the emissions of the 
two electrons in the pair is shown as a split of the wavepacket. 
Right: The transmission probabilities Ta (dashed), T!j (dot- 
ted) and Tg (solid) as a function of Tb [T v = cos 2 (<^r,)], giving 
optimal violation of the Bell inequalities for dephasing param- 
eters 7 = 1 and 7 = 0.3. 



emitted into the same normal lead, only contributes [l3| 
to the cross correlator at order T 4 and can thus be ne- 
glected. From Eq. @ we then find that S a fj((f)A, <Pb) = 
PoPaf3{4>A,4>B), where P a jj is given in Eq. (JSJ), just as 
announced. We note that it is the structure of in 12- 
space, i.e. |\&i2) , that determines the angle dependence 
of P a f3, all properties of I^ab) just gives rise to the con- 
stant Po. Moreover, the calculations show that the time 
dependent correlator (5Ii a (t)5l2f3{0)) vanishes as (t c /£) 2 
for t t c . As is pointed out above, S a f3(4>A,4>B) can 
thus be inserted in Eq. |0J and be used to violate the BI 
in Eq. ©. 

Until now, ideal conditions have been considered. One 
possible source of disturbances is dephasing. In our sys- 
tem, dephasing can quite generally be expressed in terms 
of a density matrix p = [|11)(11| + |22)(22| +7(|11)(22| + 
|ll)(22|)]/2, where the off-diagonal elements, giving rise 
to the entanglement, are suppressed by a phenomenolog- 
ical dephasing parameter < 7 < 1. The correlators 
E(4>a,4>b) in Eq. (@J then take the form 

E((j> A , 4>b) = cos(2(j) A ) cos(2c/) B ) + 7sin(2<^4) sin(2</> B ). 

(10) 

By adjusting the four angles 4>A, 4>ai and <p' B we find 
the maximal Bell parameter in Eq. @ is 



S = 2v/l + 7 2 



(11) 



which violates the BI for any 7 > 0. The optimal vio- 
lation angles, all in the first quadrant, are tan(2<^i) = 
— 7Cot(</>s), tan(20^) = 7tan(0s) and tan(0s — 
<t>' B ) = sign[cos(2<^)][(tan 2 (<fe) + l 2 )/{l 2 tan 2 (0 s ) + 
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l)] 1 / 2 , where 4>s — (\>b + 4>'b can be chosen at will. The 
corresponding transmission probabilities T v — cos 2 ((/>,,) 
are shown for 7=1 and 7 = 0.3 in Fig. [21 

The BI can thus in principle be violated for any amount 
of dephasing. However it might be difficult to produce 
beam splitters which can reach all transmission proba- 
bilities between and 1. This is not a serious prob- 
lem in the absence of dephasing, 7 = 1, a violation can 
be obtained for a large, order of unity, fraction of the 
"transmission probability space" . However, in the limit 
of strong dephasing, 7 <C 1, the set of transmission prob- 
abilities for optimal violation contain transmissions close 
to both and 1, see Fig. Expecting unity transmis- 
sion to be most complicated to reach experimentally, we 
note that by instead choosing transmission probabilities 
T A = T B = 0, T' B = 1/2 and T' A < 7, the inequality in 
Eq. J2J becomes 2\ \ + < 2. This gives a violation, 
although not maximal, for all 7 -C 1. 

Apart from dephasing there are several other effects 
such as additional scattering phases, impurity scattering 
or asymmetric tunnel barriers, which might alter the pos- 
sibility to violate the BI. It turns out that all these effects 
can be taken into account by replacing 7 — > cos((f>o) 
in Eq. (| 1 1 1) - with the important conclusion that none of 
these effects will destroy the possibility to violate the BI. 

The phase factor <f>Q is the sum of possible scattering 
phases from the beam splitters [the scattering amplitudes 
in Eq. are taken real], phases ~ kpAL due to a dif- 
ference in length, AL, between the normal arms (see Fig. 

, scattering phases from weak impurities and a possible 
phase difference between the superconductor at the two 
tunnel contacts 1 and 2. As a consequence e.g. the super- 
conducting phase (in a loop geometry) can be modulated 
to compensate for the other phases. 

The factor 7' plays the same role as dephasing. One 
possible contribution to 7' is energy dependent phases 
which oscillate rapidly on a scale of eV, suppressing the 
entangled part of the current correlator. For different 
lengths of the normal arms, there is always a phase 
~ EAL/Tivf- This phase can however be neglected 
for AL <C hvp/eV, which for eV C A is fulfilled for 
AL smaller than the superconducting coherence length 
Kvf/A. Another possibility is that, due to asymmetries 
of the tunnel barriers Ti ^ T2, the amplitude for the pro- 
cess where the pair is emitted to |11) is different from the 
process where it is emitted to |22). This gives rise to a 
statejin 12-space, (ri|ll) + r 2 |22))/ v /rf+rf). In this 
case 7' = 2rir 2 /(r 2 + T 2 ). Thus, it is in principle 
possible to violate BI for arbitrary asymmetry. 

We finally point out that the constraint on single 
mode beam splitter can not easily be relaxed. A multi- 
mode beam splitter will, in the suggested setup, probably 
have a different transparency dependence for the different 
modes. This decreases the space of accessible angles and 
will eventually make the violation of the BI impossible. 

In conclusion, we have investigated a spin-independent 



scheme to generate and detect two-particle orbital en- 
tanglement in a mesoscopic normal-superconductor sys- 
tem. The cross-correlator between the currents in the 
two leads depends in a nonlocal way on transparencies of 
beam splitters in the two leads. These nonlocal correla- 
tions arises due to the entanglement of the injected pair. 
For appropriate choices of transparencies, the correlators 
give a violation of a BI for arbitrary strong dephasing. 
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with Novel Electronic Properties. 
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